/* 
 * File:   utils.h
 * Author: jonk
 *
 * Created on 23 de noviembre de 2011, 17:28
 */

#ifndef UTILS_H
#define	UTILS_H

void histogram(int ndata, double *data, int nbins, double xini, double xfin, double *hist)
{
    int i, id;
    double step = (xfin-xini)/double(nbins);
    /*Inicializacion*/
    for(i=0;i<nbins;i++) hist[i]=0.0;
    
    for(i=0;i<ndata;i++)
    {
        id = (data[i]-xini)/step;
        hist[id]+=1.0;
    }
    for(i=0;i<nbins;i++)
    {
        /*Normalizacion*/
        hist[i]/=double(ndata);
    }
}

#endif	/* UTILS_H */

